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Abstract 

Using the recently proposed ggrf effective vertex, we investigate the pro- 
duction of rf from gluon fusion in polarized pp collisions. We show that by 
measuring All i n v' production, one can extract the polarized gluon distri- 
bution AG(x, Q 2 ) at Q 2 ~ 1 GeV 2 and in a wide range of x. 
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1 Introduction 



The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
will collide polarized protons at a center of mass energy of ^/s = 250 GeV in very 
near future (for a review of the spin program at RHIC, see [p]]). One of the most 
important goals of the experiment is to measure the poorly known polarized gluon 
distribution in a proton. In general, polarized distributions will be measured at 
RHIC in a more extended kinematic region and with higher accuracy than the 
previous experiments. 

Even though Deep Inelastic Scattering of electrons on protons is the best process 
in which to measure QCD structure functions and the parton distribution functions, 
it is also possible to extract the parton distribution functions from Drell-Yan, di- 
rect photons, J/V's, etc. Recently, we proposed [0] measuring 77' in (unpolarized) 
proton-nucleus collisions as a way to extract the gluon distribution function in nu- 
clei. This process allows one to directly (without any deconvolution) extract the 
nuclear gluon distribution function. Furthermore, the accessible kinematic region 
in x in this measurement is orders of magnitude smaller than other processes in 
hadronic collisions. One can also investigate the possibility of the restoration of 
Ua(X) symmetry at zero temperature due to high gluon density effects in nuclei 0. 
Production of rf in heavy ion collisions was considered in ||. 

In this short note, we extend our previous work on rf production in proton- 
nucleus collisions to polarized proton-proton collisions at RHIC and estimate the 
double spin asymmetry A LL f\. 

We show that the predicted asymmetry is large enough to be experimentally 
measurable in a wide range of rapidity. We predict the asymmetry All to be 
~ 0.001 at mid rapidity using the available parameterization of the polarized gluon 
distribution functions. Alternatively, by measuring the asymmetry All m the mid 
rapidity region, one can probe the polarized gluon distribution function at x ~ 0.004 
and Q 2 = 1 GeV 2 . 

2 rf production cross section 

The rf meson is much heavier than its fellow pseudo-scalar mesons, n, K and 77. 
This fact has intrigued the physics community for decades. It is now understood 
that rf gets most of its mass through instantons and breaking of Ua(1) symmetry 

1 The possibility of using rf at low P t to probe the polarization of gluons inside the proton was 
first mentioned in Ref.Q. 
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due to the triangle anomaly || [7J. Atwood and Soni used this anomaly to propose 
the following effective vertex for gluon-gluon-77' coupling in |8j 

T a f (p, g, P) = P(p 2 , q 2 , P 2 ) 6 ab e^ A7 ^ q» e x (p, a) e%, (3) (1) 

p and g are the momanta of the gluons while P is the momentum of the produced 7/. 
Here, e\(p, a) is the polarization vector of a gluon with momentum p and helicity 
a = ±1. The Kronecker delta 5 a b indicates that only the color singlet combination 
of the gluons contribute. In Ref.||, the form factor H(p 2 ,q 2 ,P 2 ) was estimated to 
be 1.8 GeV -1 in the limit where the incoming gluons and the produced rf are on- 
shell. Using this effective vertex, it is easy to calculate the differential cross section 
for gg — > rj' which is given by || 



1 ab,af3 



*ZU = n J ^flW^P - p - q)j£^r (2) 
4y/(j> ■ q) 2 [2ixy2E P 

with p 2 = q 2 = and P 2 = M 2 ,. 

We will now consider the collision of longitudinally polarized protons at RHIC. 
The double spin asymmetry A LL is defined as (for an introduction to spin effects in 
high energy collisions, see || [TO], [Til ) 



_ da ++ - da+- = dAa 
da ++ + da+- 2da 

where da ++ denotes the cross section where both protons have their spins parallel 

to their momenta while da^ denotes the cross section when one proton has its spin 

anti-parallel to its momentum. The unpolarized proton-proton differential cross 
section is denoted by da. 

The differential cross section dAa is related to the polarized gluon distribution 
function through 

dAa pp ^' x = J dx 1 dx 2 AG(x 1 ,Q 2 f )AG(x 2 ,Q 2 f )dAa"^ v ' (4) 

where x\ and x 2 are the momentum fractions of the incoming gluons, Q 2 is the fac- 
torization scale and dAa = da ++ — da + ^. Here da ++ is the differential cross section 
for scattering of two positive helicity gluons averaged over the color. Explicitly, 

\ c J a,b 
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and similarly for da + _. However, since rj' is a pseudo-scalar, it is clear that only 
positive helicity gluons will contribute. In other words, 

dAa 99 ^' = da ++ (6) 

Calculation of IT^" 1 "! 2 is straightforward and gives 

uri 2 = sj£j£ (7) 

where we defined H = H(0,0, M 2 ,). The differential cross section for inclusive rj' 
production in polarized proton-proton collisions is thenQ 

dAa ^' X = x+ AG(x + , Q)) x_AG(x_, Q)) (9) 

dx L 32 x E 

where xi — 2P z /y r s and xe = 2E P /y/s are the momentum and energy fractions of 
the produced rj' respectively and 

x E ±x L E P ±P Z 

x± " —T~ = (10) 

Using (^) in (|3|) leads to the following result for the asymmetry A LL 

_ AG(x + ,Q})AG(x.,Q}) 
LL G(x + ,Q})G(x.,Qj) [ } 

Equation (|TTJ) is our main result. It simply relates the measured double spin asym- 
metry to the polarized gluon distribution in a proton. 



3 Results 

To estimate the cross-sections and the theoretical uncertainties, we use following 



two sets of gluon distribution functions. The first set consists of the MRST99[12 



2 The unpolarized cross-section calculated in Ref. M is 
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^x + G(x+,Q 2 f )x^G(x-,Q 2 f ) (8) 



dxL 64 xe 

which is the average of the polarized cross-sections da^T** 1 X /dx^ and da^Z* 71 /dxi,. 
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Figure 1: The gluon distribution function xG(x, QV) from MRST99 (solid line) and 

n 



GRV98 (dashed line). Both are calculated at Q 2 = 1.25 GeV 2 



unpolarized gluon distribution function and LSS2001 El polarized gluon distribu 



tion function which is based on MRST99. The second set consists of GRV98[|b4 
unpolarized gluon distribution function and GRSV2000[[I5| polarized gluon distri- 
bution function which is based on GRV98. The minimum Q 2 allowed by MRST99 
is 1.25 GeV 2 and that's what we set our value of Q 2 to be. 

Using MRST99 the cross section integrated between 10 -5 < xl < 0.1 is 

CT = i(a ++ + ct+_) = 0.21 mb (12) 

while using GRV98 gives 

a = i(a++ + ct+-) = 1.0 mb (13) 

in the same range. This difference in the cross-section is mainly due to the difference 
of the two gluon distribution functions at small x. For illustration, we plot xG(x, Q 2 ) 
for both sets in Figure ([!]). 
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From the thermal model as well as exclusive reaction studies TE] one would 
expect the rf cross-section to be a few% of the pion production cross-section, or 
0(0.1 mb). However, this does not necessarily imply that MRST99 is better suited 
for our purpose. To the authors' knowledge, there is no measurement of rj' inclusive 
production cross-section. By measuring the rj' inclusive production cross-section, 
one can better constraint the behavior of the gluon distribution function at small x. 

The polarized gluon distributions from LSS2001 and GRSV2000 are also very 
different. Using LSS2001, we get 

Aa = (ct ++ - a+_) = 0.97 /xb (14) 

integrated in the 10" 5 < xl < 0.1 range while using GRSV2000 gives 

Aa = (a ++ -ff + _) = 0.17/ib (15) 

in the same range. These results indicates that the asymmetry All can range from 
1.7 x 10~ 4 to 4.6 x 10~ 3 or can differ by a factor of 30 or more depending on the 
choice of current gluon distribution functions. Therefore measuring rj cross-section 
and the asymmetry accurately can greatly help constrain different gluon distribution 
functions currently in use. 

Figure (0) shows All as a function of rf meson momentum fraction xl at y^i = 
250 GeV. Again we use the two sets of gluon distribution functions as described 
above. The solid curve in Figure (|2[) uses Set 1 (MRST99+LSS2001) and the dashed 
curve uses Set 2 (GRV98+GRSV2000). The difference between two sets is quite 
large. The values of All at small x differ by a factor of almost 40. 

One notable feature of this figure is that for xl < 0.001, All is constant. To 
understand the plateau, we note that 



fa? L + 4M%/s±x L 
x± = 1 (16) 



If xl <C M v > I a/s, then x + ~ a;_ ~ M v i / y/s and hence 

(17) 



^ ( AG{x = M„/y/8) V 



will remain almost constant. The value of the asymmetry All at y — of course is a 
function of the collision energy. The input parametrization used for the unpolarized 
gluon distribution function for both LSS2001 and GRSV2000 is of the form 

xAG(x) =Cx a (l-x) b xG{x) (18) 
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Figure 2: The asymmetry A LL as a function of rj' Xl calculated at ^/s = 250 GeV 
and Q) = 1.25 GeV 2 . The upper curve is calculated using MRST99 and LSS2001 
parton distribution functions. The lower curve is calculated using GRV98 and 
GRSV2000 sets. 

If our Q 2 j = 1.25 GeV 2 is not very much different from the input Ql of the distribu- 
tion functions, the asymmetry should be given by 



A L l(xl) 



C 2 x a + (l-x + ) b x a 4l-x_) 



b 




(19) 



where we used x + X- = M 2 , 



/s. Near xl = 0, both x + and x_ are small so that 




(20) 



Therefore measuring the asymmetry at y = for various energies can tell us what 
the exponent a should be without knowing the detailed form of xG(x). For LSS2001, 



this exponent turned out to be about 0.6 and for GRSV2000, it is about 0.8. The 
reason that Set 2 curve deviates from a straight line may be explained by the fact 
that the input Ql of set 2 (0.4 GeV 2 ) is not so close to our Q 2 f = 1.25 GeV 2 . 
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Figure 3: The asymmetry A LL at the mid-rapidity calculated at Q"j = 1.25 GeV 2 . 
MRST99 and LSS2001 are used for the upper curve and GRV98 and GRSV2000 
are used for the lower curve. The solid lines are linear fits. For MRST99+LSS2001, 
the slope turned out to be —1.27 corresponding to a = 0.633 in Eq. (|18D . For 
GRV94+GRSV, the slope is -1.65 corresponding to a = 0.827. 

Even though the detectors at RHIC can only go up to the maximum x of about 
0.01, it is worth noting the behavior of An for larger values of x in Figure (0). 
When 1^ > M v '/^/s, x± can be approximated to be| 

M l> , K , , 

x + « xl H and a;_ w — - (21) 

SXl SXl 
3 For RHIC at = 250 GeV, M v >/^/s = 0.0038. For LHC at y/s = 5500 GeV, M^/y/s = 
0.00017. 
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and the asymmetry becomes 



In LSS2001, the exponent 6 is argued to be 0. In other parametrizations, b is 
substantially different from 0. For instance in GRSV2000, b ranges from 4 to 7. 
Therefore even at moderately large xl, the behavior of All as a function of xl 
depends strongly on the exponent b. In this way, very forward 7/'s can distinguish 
between different parameterizations of the unpolarized gluon distribution function. 



4 Conclusion 

In this paper, we showed the dependence of A L l on rf xl as well as the (pp) center 
of mass energy of the collison using ([11]). Measuring All i n this way directly probes 
the gluon distribution function in the proton. Had we known the unpolarized gluon 
distribution function well, A L l could be used to map out the polarized gluon distri- 
bution. As it happens, currently the gluon distribution is not very well constrained 
by the past experiments in the small x region. Measuring rj' at high energy can 
greatly constrain the gluon distribution in this region. 

There are several caveats to our results. First and perhaps the most important 
one theoretically, is that M v > ~ lGeV is not very large and therefore higher order 
(in a s ) corrections can be quite sizable. This is somewhat "cured" by introducing 
a K factor which is supposed to take higher order corrections into effect. In our 
case, since we are calculating the ratio of cross sections, these effects should largely 
cancel out. However, without a real calculation of higher order effects, there is no 
way to know for sure. Another complication arises due to effects such as initial 
state radiation of gluons, etc. which would lead to the produced particle having 
a transverse momentum. These effects can be resummed for few limited processes 
but are mostly included phenomenologically by introducing an intrinsic k t for the 



incoming partons fll7| . We intend to do this in the near future and investigate the 
dependence of All on the transverse momentum of rj'. However, if we consider p t 
smaller than the intrinsic k t scale, our analysis should still apply. 

Another potential problem is that we have neglected possible contribution of an- 
nihilation of quarks and anti quarks into rj's. We believe we can avoid this problem 
here for two reasons. First, the values of x considered are quite small where glu- 
ons are the dominant partons. Therefore contribution of sea quarks to this process 
should be of order of a few percent which we can safely neglect. Furthermore, if 
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the qqrj' vertex is of the derivative type[[nj, then in the on-shell limit the quark 
anti-quark annihilation matrix element is 0{m q /M v t) smaller than the gluon fu- 
sion matrix element. Hence, in the small x range the gluon fusion process should 
dominate over quark anti-quark annihilation process. 

In summary, we have investigated the production of rj' in polarized pp collisions 
and estimated the double spin asymmetry variable All using the available polarized 
gluon distribution functions. Alternatively, measuring All at polarized proton- 
proton collisions at RHIC will enable us to directly extract the polarized gluon 
distribution function at Q 2 ~ 1 and different values of x. 
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